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Abstract 



Due to considerable recent interest in the use of density matrices for a wide va- 
riety of purposes, including quantum computation, we present a general method 
^ \ for their parameterizations in terms of Euler angles. We assert that this is of more 

' fundamental importance than (as several people have remarked to us) "just an- 

!^ , other parameterization of the density matrix." There are several uses to which this 

'sj- \ methodology can be put. One that has received particular attention is in the con- 

■ struction of certain distinguished (Bures) measures on the (n^ — l)-dimensional 

^ . convex sets of n x n density matrices. 

^: 
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1 Introduction 



The density matrix [1] — having origins in early (independent) work of Lan- 
dau and von Neumann — has proved to be a very useful concept in physics. 
Researchers have devoted substantial efforts in describing the spaces defined 
by density matrices [2], in using them to analyze the separability of quantum 
systems [3,4], in comparing information-theoretic properties of various proba- 
bility distributions over them [5], as well as studying the question of parallel 
transport in this context [6]. There is, of course, a quite straightforward man- 
ner in which to parameterize density matrices — simply in terms of the real 
and complex parts of the entries of these Hermitian matrices. A more indirect 
approach, but one of particular interest, relies upon Euler angle parameteri- 
zations of the special unitary matrices. It has been known for quite some time 
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that one can parameterize density matrices with the help of such angles. To 
do so, we take a diagonal density matrix, p, which represents our quantum 
system in a particular basis. We then perform a unitary transformation in 
Hilbert space that takes p to an arbitrary basis. This unitary transformation 
can be assumed to have a determinant equal to one since an overall phase 
does not affect the physics. Now when we apply the unitary transformation, 
we operate in the following manner, 

p'^UpU\ 

where U G SU (n) for an n-state system. The diagonal matrix p can be pa- 
rameterized by its n — 1 eigenvalues, and the unitary transformation by — 1 
variables. At this point, however, we have parameterized p' using too many, 
that is n^+n — 2 parameters, since the density matrices for the n-state systems 
comprise only an (n^ — l)-dimensional convex set. A primary objective here 
will be to demonstrate how such an "over-parameterization" can be avoided, 
by the elimination (from the unitary transformation) of n — 1 of the set of 
+ n — 2 parameters. 

Below, we will show that the density matrices for the 2-statc and 3-state sys- 
tems can be conveniently and insightfully parameterized in terms of Euler 
angle coordinates. In doing so, we take full advantage of the group properties 
associated with the unitary transformation. In particular, we obtain useful 
(Bures) measures on the spaces of density matrices, which can be used in 
the integration of functions over spaces and subspaces of quantum systems. 
Specifically, in Section 2.1 we exhibit the Euler angle-based parameterization 
of the 2-state quantum systems. In Section 2.2 we present the Euler angle- 
based parameterization of the 3-state systems. In Section 2.3 we discuss their 
straightforward generalization to n-state systems. In Section 3.1 we show how 
to use this to write down the Bures measures on the space of density matri- 
ces for 2-state systems and in Section 3.2 its counterpart for 3-state systems. 
Finally, we discuss ongoing work which is intended to yield the exphcit gen- 
erahzation to 4-state systems. This, we anticipate will be highly useful in 
analyzing the (fifteen-dimensional) space occupied by pairs of qubits (cf. [7]). 



2 Pcirametrization of Density Matrices 

2.1 2-state Density Matrices 

In this section we review certain properties of density matrices for the 2- 
state systems and demonstrate how to parameterize the 2-state density matrix 
in Euler angles. This should make the extension to 3-state systems in the 
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following section more transparent. 



Much is known about 2-state density matrices. One of the more obvious prop- 
erties is that a 2-state system can be spanned by two pure state density ma- 
trices of the following form: 

fio\ (oo\ 

Pi = ' ^2 = . (1) 

(These can be denoted by antipodal points on a two-sphere, making use of 
the elegant Riemann sphere representation [11].) For a mixed state diagonal 
density matrix, we write the linear combination of these two density matrices 
as 

P = EPi«'' (2) 

where c^* = 1- The a* may be represented in several ways. One common way 
is to let = a and = 1 — a. What we wish to emphasize in this paper is 
that they are better represented by — cos^ 9 and = sin^ 9. Why this is 
true will be more obvious in the next section. 

To take p to an arbitrary configuration (basis) we act with a unitary transfor- 
mation U G SU{2), in the manner 

p' ^ UpU-^ = UpUl (3) 

One may note that an n-state density matrix should have — 1 parameters, 
whereas this one appears to have 4, that is 3 (from SU (2))-|-l (from the diagonal) . 
However, in the Euler angle parameterization of SU{2), given by 

where the as are the Pauli matrices, we see that the parameter 7 drops out 
since (T3 is diagonal and so the matrix exponential e'^'^^'^ commutes with p 
leaving precisely three parameters, as is appropriate and desired. 

An important part of this parameterization is the ranges of the angles. For 
the 2-state case these are well-known (up to normalization given in eg., [12]). 
They are 

0<Q;<7r, 0</3< 7r/2, < ^ < 7r/4. (5) 
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2.2 3-state Density Matrices 



For 3-state density matrices we can achieve a similar parameterization. The 
pure states are spanned by 



Pi 



^10 0^ 


yO Oy 



P2 



'^O 0^ 
1 

yO Oy 



P3 



^0 0^ 




(6) 



We can then obtain a mixed state by the hnear combination J2PiO'\ with 
X] a* = 1. Here again we use the squared components of a sphere. In this case, 
it is a two-sphere, thus for a 3-state system we take a generic mixed state to 
be 



P 



^ cos2 01 sin^ 02 

sin^ Oi sin^ 62 

cos^ 02 



\ 



(7) 



We then can take this to an arbitrary configuration (basis) by the process 



P^p' = UpU\ 



(8) 



where U e SU{Z). In the Euler angle parameterization of 5^/(3), U is given 
by [13] 



JJ — giA3agiA2/3gjA37giA50giA3agjA2bgiA3CgiA8</>/A/3 



(9) 
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where the Gell-Mann matrices have been used; 



A4 



'^0 1 0^ 



1 




^0 1^ 





1 



'^O ^ 



A7 = 



-i 

i 



A, 



^0-^0^ 



i 
00 
-i 



1 
1 
1 
0-2 



A, 



1 
0-10 




^0 0^ 



1 
1 



(10) 



Since A3, and Ag appear on the right in (9), they commute with the diagonal 
matrix p and drop out of the parameterization. Thus we are left with 3^ — 1 = 8 
parameters for the 3-state case, as we should expect [2]. 

As in the 2-state case, the ranges of the angles in the parameterization are 
very important. For the case of three states we have 



0<Q;,7,a<7r, < P,9,b < 7r/2, 
0<0i< 7r/4, and < ^2 < cos-^(l/V3). 



(11) 
(12) 



Here one should note that using a, b and 1 — a — 6 for the diagonal elements 
of the density matrix would give rise to a domain of integration that is non- 
rectangular. For density matrices of higher dimensional systems this could be 
extremely akward. Here however, we see that the domain is indeed a rectan- 
gular solid with the final angle for a density matrix of n-dimensions having 
the range < 9n-i < cos~^ {1 / ^/r^) . 



2.3 n-state Density Matrices 



Though there is now no explicit parameterization of SU{n) in terms of Eu- 
ler angles, the generalization is rather obvious. The difficulty in manipulat- 
ing the one-forms to calculate the volume elements increases as so un- 
less some substantial progress is made in getting computers to manipulate 
trigonometric functions symbolically, it appears that progress will be quite 
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slow. It is anticipated that the density-matrix-parameterization program ad- 
vanced here will not be achieveablc for groups much higher in dimensionality 
than dim{SU{4:)) = 15. Perhaps, however, with sufficient interest and effort it 
could be accomphshed for the Lie group SU (8), in order to effectively describe 
the 63-dimensional space of three entangled qubits. 



3 Bures Measures on the Convex Sets of Density Matrices 



3.1 2-state Systems 



The parameterization of the space of density matrices by Euler angles imme- 
diately leads to a particularly simple procedure for integrating over the space. 
A natural (Bures) measure on the space is given by the product of the (Hall) 
measure [15, eqs. (24, 25)] [16] on the space of eigenvalues (Ai, . . . , A„), that 
is 

<i«=^^^n4%^ (13) 

(Ai . . . An) 2 ^3 + 



and the measure on the the space of unitary matrices, the (modified/truncated) 
Haar measure. The geometry of the three-dimensional space is obvious from 
the parameterization. We have a solid section of the two-sphere and the space 
SU{2)/U{1) = 5"^. In other words, our measure is given by 

DV^dux dS^, (14) 



where DV is the Bures measure on the space of density matrices, dS"^ is the 
standard measure on S^, ie. dS"^ = d[G/H), with G — SU{2) and H — U{1) 
[8] and du is given by (13) with n = 2. 



3.2 3-state Systems 



The measure on the eight-dimensional space of density matrices describing 
3-states is a generahzation of the previous case. We have 

DV ^duxd{G/H), (15) 



where DV is the Bures measure on the space of density matrices, with G ~ 
SU{3) and H = U{1) x U{1) [8] and du is given by (13) with n = 3. The 
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ranges of the Euler angles are given in Section 2.2. This is a new result. In 

the past the ranges of the angles were not properly specified when using the 
Haar measure. Additionally, the problem of overparameterization [3] was not 
clarified until recently [10]. 

In [16], it has been shown how the Bures measures for the 2- and 3-state 
systems can be normalized to form probability distributions. 



4 Applications / Conclusions 

The point of this note has been to clarify that: (1) the Euler angle parameter- 
ization eliminates any naive over-parameterization ([3]) of the density matrix 
using a unitary group representation; and (2) doing so helps in providing a 
natural (Bures) measure on the space of density matrices, which incorporates 
the (truncated) Haar measure on the group manifold of SU (n) . 

So this is not "just another parameterization of the density matrix", while 
proving very useful for many calculations. We find that people often do not 
appreciate the point of view of Feynman, who wrote that "every theoretical 
physicist who is any good knows six or seven different theoretical representa- 
tions for exactly the same physics" [9]. 

In addition this gives rise to a rather natural (Bures) measure on the space 
of density matrices, using which one can integrate functions of quantum- 
mechanical interest (such as the von Neumann entropy, — Trp log p) over spaces 
and subspaces of the group with relative ease. Much of the geometry of the 
space can then be seen in terms of the moduli of the group space [10], a point 
apparently not appreciated by a number of readers of [10]. 
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